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ABSTRACT 

Observations of white dwarfs in the globular clusters NGC 6397 and Omega Centauri 
indicate that these stars may get a velocity kick during their time as giants. If the mass 
loss while on the asymptotic giant branch is slightly asymmetric, the resulting white dwarf 
could be born with such a velocity kick. These energetic white dwarfs will impart their 
excess energy on other stars as they travel through the cluster. A Monte-Carlo simulation 
of the white-dwarfs kicks combined with estimate of the phase-space diffusion of the white 
dwarfs reveals that as the white dwarfs equilibrate, they lose most of their energy in the 
central region of the cluster. They could possibly augment the effect of binaries, delaying 
core collapse or increasing the size of globular cluster cores. 
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1 INTRODUCTION 

Several phenomena give circumstantial evidence that white 
dwarfs may receive a velocity boost (or kick) at birth. The 
observed rotation rates of white dwarfs could result from mild 
kicks (~ km/s) generated by asymmetric and off-centered 
( winds toward t he end of th eir time on the asymptotic giant 
branch (AGB) (|Spruitll998l h It is unclear at present how white 
dwarfs could actually obtain such a kick; however, given the 
evidence it is natural to explore the consequences of white- 
dwarf kicks. Mild kicks m ay explain the possible l ack of white 
dwar f s in open clusters ( Fellhauer et al.l 20031; Weidemarml 
1 19771 ; iKalirai et al.l lioOll ). Most directly, iDavis et all (|2008l ) 
observed that the young white dwarfs in NGC 6397 had a more 
expansive radial dis tribution than e i ther t heir progenitors or 
older white dwarfs. ICalamida et alj (|200&f ) found similar but 
weaker hints in Omega Centauri. 

Without a kick young white dwarfs would have a velocity 
distribution nearly equal to that of their more massive progen- 
itors on the main sequence. In this case the kinetic energy of 
these white dwarfs is much less than equipartition; therefore, 
as their velocity distribution relaxes they cool the rest of the 
cluster. If on the other hand white dwarfs r eceive a substantia l 
kick at birth as observations may indicate (|Davis et alj|2008l h 



young white dwarfs may heat the rest of the stars in the clus- 
ter. Neutron stars that form in the cluster may get a very 
large velocity kick as well; however, two factors make their 
contribution less important. The first is that neutron stars are 
relatively rare, and the second is that neutron stars travel so 
quickly that most escape the cluster before imp arting much 
energy to the rest of the st ars ([Pfahl et al.ll2002f ). In contrast 
with the neutron-star result IHevll ( 20071 ) fo und that mos t of the 
kicked white dwarfs remain in the cluster. IHevll (|2008bh found 
that the energy input from the white-dwarf kicks around the 
half-light radius is about half that from binaries at the present 
day and dominated in the first half of the life of the cluster. 
This letter examines where in the cluster white dwarfs dump 
their excess energy from the kicks. 



2 CALCULATIONS 

A globular cluster is often model led wit h a lowered i s other- 
mal profile (or King m odel) (|Michiel 1 19631 ; iKind Il966l ; 
iBinnev fc Tremaindll987lh 
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where 



and is the gravitational potential, a is 



the velocity dispersion of the cluster stars, and pi is a char- 
acteristic number density. The distribution function depends 
only on the energy, —me, a constant of the motion; therefore, 
it is constant in time as well. 

As the stars interact with each other, their kinetic energy 
approaches equiparti tion such tha t m^cr 2 = m^-cr 2 for masses 
mi not equal to rrij (Spitzer 1987). The most massive main- 
sequence stars in the cluster are the progenitors of the white 
dwarfs so so they will typically have <tto < "cluster, where 
^duster is the mean velocity dispersion of the cluster. To model 
the cluster, the turn-off star s and the wh ite-dwarf kicks, we use 
the best-fitting model from iHevl (|2008al ) . Specifically, we take 
""cluster = 1 for the cluster stars, <jto = 0.5 for the turn-off 
stars, the central potential to ^(O) = 6 and the mass of the 
white dwarfs to be 1.5 times the mass of the typical cluster 
member. We model the kick as a Gaussian in velocity with a 
dispersion of <Jk = 0.92. This means that the kick is similar 
to the velocity dispersion of the cluster. Although this is rea- 
sonable for clusters such as NGC 6397 with <7 c iustcr ~ 5 km/s, 
it is more difficult to achieve for clusters such as 47 Tuc and 
Omega Cen with <7 c i us tcr > 20 km/s. 
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2.1 Relaxation from a kick 



The rate that a particular st ar gains energy as it pa sses 
through the cluster is given by (|Binnev fc Tremainelll987l ) 
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where m and v are the mass and velocity of the star, m a and v a 
are the mass and velocity of the other cluster members, In A is 
the Coulomb logarithm and f a {v a ) is the phase-space density. 
The first term is the energy that a particular star gains from 
encounters from faster moving stars. For a King model it is 



V a fa(Va)dv a = pi (27TCr 2 ) (J 2 fe' /<T - lj 
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where e = \& — |u 2 for the particular star, a is the dispersion 
for the cluster as a whole. The second term is the energy that 
the star loses from encounters with slower moving stars, 



Figure 1. The power imparted to the cluster from young white 
dwarfs within a given radius as a function of the enclosed cluster 
mass. The upper curve traces the value for white dwarfs with a kick. 



3 RESULTS 

To explore the distribution of the energy from the kicks to 
the rest of the cluster members, we created a Monte-Carlo 
realisation of 100,00 newly bo rn white dwarfs both with and 
without kicks (see iHevll [2007, for further details) . For each 
realisation the white dwarfs are sorted by radius and the total 
energy deposited into the cluster stars is summed by radius 
resulting in the cumulative distribution of power depicted in 
Fig. [1] The upper curve traces the result for white dwarfs 
with kicks. The white dwarfs are travelling too fast for their 
mass so they deposit the excess energy into the cluster. The 
lower curve shows the result for white dwarfs without a kick 
— here, the white dwarfs are born with the velocity dispersion 
of their progenitors who are significantly more massive than 
themselves; therefore, the white dwarfs sap energy from the 
cluster stars as they heat up. The exact locations of the two 
curves depend on the assumed value of the ratio between the 
white-dwarf mass and that of the cluster stars (the figure gives 
the result for 1.5; as the ratio increases both curves move up). 
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Because we know the position and velocity of each of the white 
dwarfs in the simulation, it is straightforward to calculate the 
energy loss and gain from these equations. 



The distance between the cumulative energy deposition 
rates is the net effect of the white-dwarf kicks. This net effect is 
rather robust because regardless of the assumed ratio between 
the white dwarf masses and the rest of the stars, white dwarfs 
are typically created near the centre of the cluster because 
their progenitors are massive. Furthermore, even with a kick, 
the white dwarfs will return to the radius of their birth until 
their orbits relax. At this small radius, the speed of the white 
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Figure 2. Net increase in energy deposition from the kicks within 
a given radius as a function of radius and the mass enclosed with 
that radius. The dashed curve panel traces the cumulative mass 
distribution. Both the cumulative (solid curve) and the differential 
power distributions (dot-dashed curve) are shown. The tidal radius 
of the cluster is r< = 0.77, the core radius is r c = 0.04, and the half- 
light radius is r h = 0.11. Ninety percent of the energy is deposited 
within a core radius of the centre of the cluster. Beyond a radius 
of 0.14, the kicked white dwarfs begin to lose energy to the cluster 
stars. 



dwarf is largest and so is the density of the rest of the stars; 
therefore, at small radii the white dwarfs lose most of their 
excess energy. Fig. [2] demonstrates this fact. For the assumed 
model about ninety percent of the energy of the white-dwarf 
kicks is deposited within the core radius, and the white dwarfs 
actually remove energy from the outer regions of the cluster. 

During the collapse of a globular cluster core, two-body 
interactions result in the transfer of energy from the dense cen- 
tral region of the cluster to the outer regions; the core shrinks, 
and the envelope of the cluster expands. The white dwarfs 
deposit their energy in precisely the opposite way, supplying 
energy to the core at the expense of the outer regions; there- 
fore, white-dwarf kicks may retard core collapse or increase 
the size of cluster cores. 



4 CONCLUSIONS 

After a globular cluster forms, the negative heat capacity 
of the cluster's self-gravity inexorably draws the cluster to- 
ward core collapse, unless another energy source is present. 
The central region or core will collapse to very large den- 



sities within a few relaxation times while the outer regions 
expand, and stars evaporate from the system. The standard 
energy s ource that is invoked to retard the collapse are bina- 
ries (e.g iGoodman fc HutJ [19891 ). After the core has already 
reached high densities, binaries necessarily form through tidal 
and three-body interactions preven ting the core from forming 
a singularity (e.g lHurlev et al.ll2007T ). 

iV-body models with a primordial binary fraction of a few 
percent can avoid the complete collapse of the core but the core 
does contract dramatically over the age of the cluster. If all 
old globular clusters had compact cores, this scenario would 
not pose a problem; however, several old, presumably relaxed 
clusters have large cores wit hout an obvious energ y source 
to explain the observed state (|De Marchi et af1l2007l ). Several 
possibilities have been invoked to explain this: a very high ini- 
tial bi nary fraction approaching unity (however, iHeggie et al.l 
|2006l ) found that the effect of binaries saturates at about ten 
perce nt), an intermediate-m ass black hole (e.g. iTrenti et al.l 
120071 ) or a black-hole binary (iMackev et al.ll2008l ) at the cen- 
tre of the cluster. As for the white-dwarf kicks, the evidence 
for these possibilities is circumstantial. 

There are several lines of evidence that white dwarfs get 
a mild velocity kick at birth. Furthermore, the total power de- 
posited into the cluster by the w hite-dwarf ki cks may be com- 
parable or larger than binaries l|Hevll l2008bl : although these 
calculations should be generalised to the core for complete- 
ness) and the energy is deposited in such a way that it directly 
opposes core collapse, so these kicks may provide a explanation 
for the fraction of globular cluster that have large cores today. 
We do however see a variety of core sizes including collapsed 
cores, so the combination of white-dwarf kicks and binaries 
must not always be effective or perhaps the clusters are at 
various stages of their evolution. Explaining these differences 
in the white-dwarf kick model is difficult as one would expect 
that all globular clusters should exhibit the effects of the kicks. 
On the other hand the black hole models can easily explain 
these differences; some clusters could have a large black hole 
at their centres and others lack one. 

The Monte-Carlo model present here is illustrative, but 
rather rudimentary. Two obvious avenues for further study 
are the inclusion of w hite-dwarf kicks in iV-body models 
(e.g. iHurlev et all I2007T I and Monte Carlo simulations (e.g. 
IHeggie fc Gierszll2008l ) — either of these avenues would allow 
us to include both binaries and white dwarf kicks together con- 
sistently to get a more robust estimate of their relative impor- 
tance. Furthermore, one could use a more realistic multi-mass 
King model for the background stars in the cluster to perform 
an analysis similar to that presented here. One would expect 
that that results could differ quantitatively but the qualita- 
tive result the white-dwarfs kicks dump energy in the central 
regions of the cluster and that this energy may be comparable 
or larger than that from binaries is robust. 
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